Electron spin relaxation in n-type InAs quantum wires 
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We investigate the electron spin relaxation of n-type InAs quantum wires by numerically solving 
the fully microscopic kinetic spin Bloch equations with the relevant scattering explicitly included. 
We find that the quantum-wire size and the growth direction influence the spin relaxation time by 
modulating the spin-orbit coupling. Due to inter-subband scattering in connection with the spin- 
orbit interaction, spin-relaxation in quantum wires can show different characteristics from those in 
bulk or quantum wells and can be effectively manipulated by various means. 

PACS numbers: 72.25.Rb, 73.21.Hb, 71.10.-w 
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I. INTRODUCTION 

Spintronics continues to attract interest because of 
potential applications to information technology, but 
also because it has greatly improved our understand- 
ing of the role of different spin-dependent interaction 
mechanisms, i.e. of the spin-orbit coupling, for many- 
electron systems jii^ Several spintronic devices have been 
proposed that manipulate the carrier spin via spin- 
orbit coupling (S0C)4i^i^ In recent years, progress 
in nanofabrication and growth techniques has made 
it possible to produce high-quality quantum wires 
(QWRs) and investigate physics in these semiconduc- 
tor nanostructuresi^'^'^ dOdid^ The energy spectrum 
of QWR systems with strong SOC has been stud- 
ied experimentallyS2,i^ and theoretically-i^ii^ii^iiiii^ii^ 
For p-type QWRs, most of the works concentrate 
on the energy spectrumjiii^ but for n-type quan- 
tum structures, many investigations have been per- 
formed with the aim of understanding the electron spin 
relaxationffi>^.25^Si28-29,30,3i 

For QWRs, the spin relaxation time (SRT) was 
measured)^!^ and calculated in the framework of 
a single-particle model^ii^i^^i^ and Monte-Carlo 
simulations Recently, Liu et al?^ investigated 
the SRT for QWRs in the (110) crystal direction, but 
systematic studies of spin relaxation in ri-type QWRs are 
scarce. A general way to obtain the nonequilibrium spin 
dynamics and the spin relaxation time in semiconductors 
heterostructures has been developed by Wu et q/. s^, 33^34 
In this approach, the momentum and spin dependent 
distribution functions are calculated using microscopic 
kinetic spin Bloch equations (KSBEs), which include 
the momentum-dependent Dresselhaus and/or Rashba 
SOC together with the effect of the relevant scattering 
mechanisms i^i^^i^ Cheng et al. applied this approach 
to study electron spin relaxation in QWR systems 
However, the influence of higher subbands was not 
included in this work. Subsequently, the effects of higher 
subbands and their coupling via Coulomb scattering 



were shown to important for the spin-relaxation of holes 
in QWRs'^^" and for the case of electrons in quantum 
wells 1^ 

In this paper, we study the influence of higher sub- 
bands on the spin relaxation in n-type InAs QWRs. This 
allows us to investigate QWRs with a wide range of 
sizes. Especially for larger diameters of the wires, the 
spin precession and the spin relaxation are expected to 
show different characteristics from narrow ones, as elec- 
trons populate more than one subband and therefore 
experience different SOCs and undergo intra-subbband 
and inter-subband Coulomb scattering. We find that the 
inter-subband Coulomb scattering can make an impor- 
tant contribution to the spin relaxation. We also study 
the influence of the growth direction of the QWRs on 
the spin relaxation. In particular, we consider QWRs 
with (001), (110), and (111) growth directions, and show 
that the SRT depends sensitively on the growth direc- 
tion, quantum-wire size and the direction of the initial 
spin polariazation. 

This paper is organized as follows: In Sec. II we de- 
scribe our model and the KSBEs. Our numerical results 
are presented in Sec. III. We conclude in Sec. IV. 



II. MODEL AND DYNAMICAL EQUATIONS 

We model the InAs QWR by a rectangular confinement 
potential, i.e., we assume infinitely high barriers at x = 
-iiax and y = ±ay, and no confinement on the z direction. 
The Hamiltonian, which describes the electronic single- 
particle states in the QWR is then taken to include the 
confinement potential Vc, the Rashba term Hn and the 
Dresselhaus term Hjj 
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Hr + Hd + Vcir) 



(1) 



The Rashba and Dresselhaus terms are the two contri- 
butions to the internal fc-dependent effective field, which 
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leads to the Dyakonov-Perel' spin dephasing mechanism. 
The Rashba term 

HR{k) = -ilf^a- • k X f 

~t~ 0'z{kxEy kySx)] (2) 

is due to the inversion asymmetry of the crystal structure. 
The Dresselhaus term is different for different growth di- 
rections. For a (100) InAs QWR, the x, y and z axes 
correspond to the [100], [010] and [001] crystallographic 
directions, respectively, and the Dresselhaus term can be 
written aSfSS, 



Hh"' = blf^{aAkx(kl-kl)]+ay[ky{kl~kl)] 



+ azikzikl - kl)]} 



(3) 



For a (110) QWR, the x, y and z directions correspond 
to the [110], [001] and [110] crystallographic directions, 
and we have 
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blf'Wx[-Mkl -kl + 2kl)] + 2aykxkykz 



+ <Jz[Mkl-kl-2kl)]} 



(4) 



For a (111) QWR, the x, y and z directions correspond 
to the [112], [110], and [111] crystallographic directions, 
and we have 



Hjj — b^-^ \^x[ "^kxkykz ~k. 
V3 
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As input for the KSBEs we use the basis {(/)„^„ ka} of 
single-particle states, which are obtained from the eigen- 
functions of + Vc{r). For a hard- wall confinement 
potential which constricts the electrons in the x and y 
directions on mesoscopic length scales, we can employ 
the envelope function approximation!^ Thus we write 
the single-particle states in the form 



(6) 



with 



i'n^,nyix,y) 



y/axa 
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sin(^)sin(^^) , (7) 



where x<y denotes the basis vectors in spin space, i.e., 
eigenstates of ctz. In the envelope function approxi- 
mation, the effective Hamiltonian acting on the single- 
particle states is obtained from Eqs. (I2][5]) by the re- 
placements kx {'>Pn^\kx\i'n') = {kx)n^,n' and kl 



{kx)n^,n'^ where kx = —id/dx. A similar replacement is 
done for ky. For our choice of the confinement potential, 
we have in particular 



(kx) = 



Aihrixrix 



axiin'x)^ - {uxY 



r(l - ^n^.n') 



[kx) = 2 — 



(8) 



with corresponding results for (ky) and (ky). From 
Eqs. ([3]) and ([8]) one can see that when and ay are 
sufficiently small and only the lowest subband in QWR 
is important, only the third term of Eq. ([3]) is not zero, 
which means the effective magnetic field contributed by 
the Dresselhaus term is along the z-direction. On the 
other hand, the third term in Eq. ^ is zero but the 
terms which proportional to ax and ay are not zero. This 
means that the effective magnetic field contributed by the 
Rashba term is in the x-y plane when the wire width is 
sufficiently small. Nevertheless for larger width of the 
wire, when higher subbands are needed, all the terms 
in Eqs. ([2|) and ^ contribute. Moreover, one can see 
that when the confinement in the x and y directions is 
symmetrical, the third term of Eq. ([3]) is zero, so that 
the effective magnetic field contributed by the Dressel- 
haus term does not contain any component along the 
z-direction. As we show in detail below, these differences 
in the spin-dependent single-particle states lead to signif- 
icantly different behavior of the SRT. 

The complete dynamical information about spin- 
dependent single-particle properties is contained in the 
spin-density matrix p. Its matrix elements are, in gen- 
eral, defined with respect to the complete set of quantum 
numbers n^,, Uy, fc, cr, but because of the translation sym- 
metry in z direction, p is diagonal in k, i.e., p = pk,ss' 
with s =^ {ux, ny,a). 

We construct the KSBEs by the non-equilibrium Green 
function method as follows:^ -■^'^ 



dp dp 

at ~dt 



dp 

coh dt 



The coherent terms can be written as 



dpk 
dt 



coh 



-i[^yQlQPk-qI-Q, Pk 
Q 



(9) 



H,{k),pk (10) 



where [A, B] ~ AB — BA denotes the commutator, and 
Q = {qx,qy,q)- Iq is a matrix in (s, s') space and can be 
considered as a form factor. Its definition reads 



Scrua2F{mi,m2,qy, ay)F{ni,n2, qx, a^OH) 



where 
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F{mi,m2,q, a) — 2mg[e*"'' cos 7r(TOi — TO2) — 1] 



1 



1 



I 



(12) 



The first term in Eq. (jTU]) is the Coulomb Hartree-Fock 
term, and the second term is the contribution from the 
single-particle Hamiltonian, i.e., Eq. ([1]) in (s,s') space, 
which includes the spin-orbit coupling terms. For small 
spin polarization, the contribution from the Hartree-Fock 
term in the coherent term is negligible^i^ and the coher- 



ent spin dynamics is essentially due to the spin precession 
around the effective internal fields described by Eqs. ([2][5]) 
The scattering contributions to the dynamic equa- 
tion of the spin-density matrix include scatterings be- 
tween electrons and nonmagnetic impurities, electrons 
and phonon, and electron and electron scatterings: 



dpk 
dt 
dpk 
dt 

dpk 
dt 



dpk 
dt 
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ph 


dt 
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J2 \UQ\^^iEs,,k-q ~ E,,,k)Ici\(X - Pk-,)TsJ-ClTs,Pk - Pk-gTsJ-QT,,il - pk)] + h.c. , 

Q,Sl,S2 



Q,si,S2,A 

- NQ^xPk-qTsJ-qTs^il - Pk)] + 5{Es^^k-q - Es^^k 

- (iVQ,A + l)pk-qTsJ-ciTs,{l - Pk)]} + h.c. , 



Es2,k + t^Q,A)[(^Q,A + 1)(1 - Pk^q)TsJ^qTs.,Pk 

^Q,a)[^Q,a(1 - Pk-q)TsJ-QTs2Pk 



'^E E V^6iEs,,k-q~Es 



Es4,k'-q)lQ 



^S2 ,k 
Q,fc' Si,S2,S3,S4 

X {(1 - P/c-g)7'si^-Ql's2PfcTr[(l - Pk')TsjQTs4Pk'-qI-Q] 

- PW-qTsJ-QTs^il - Pk)^'^[Pk'TsjQTs4{l - Pk'-q)I-Cl]} + h.C. 



(13) 
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in which Tsi^s,s' ~ <5si,s<5si,s'- The statically screened 
Coulomb potential in the random-phase approximation 
(RPA) reads'^" 



(14) 



with the bare Coulomb potential vq — Aire^ /Q^ and 

Ik+q fk 
^k+q - fife 



«(g) = i- E^qI^qI'E' 



(15) 



qm^qy 



In Eq. (fT3|) . Ni is the density of impurities, and |C/qP is 
the impurity potential. Further, |Mq^Ap and Nq,\ — 
[eicp{ujQ,\/ ksT) — are the matrix element of the 
electron-phonon interaction and the Bose distribution 
function, respectively. The phonon energy spectrum for 
phonon mode A and wavevector Q is denoted by (^q,\. 
For the electron-phonon scattering, we include the elec- 
tron LO-phonon and electron AC-phonon scattering, for 
which the explicit expressions can be found in Refs. [33I 
andli^. Note that, due to the weakness of the SOC, the 
energy dispersions Eg^k in the scattering terms, Eq. (fT3|) . 



do not include the SOC, as in Ref. [S^- This is different 
from our previous work in p-type QWR systems)^ For 
holes, the SOC is stronger because holes are derived from 
p-orbitals, which experience the spin-orbit interaction di- 
rectly, and it is therefore necessary to include the SOC 
in the single-particle energy dispersions , 35i36 ,43 



III. NUMERICAL RESULTS 

We numerically solve the KSBEs for the spin density 
matrix p, from which we obtain the dynamics of the av- 
erage spin for electrons with momentum k via 

{S)k{t) = E Pk,7i:,,ny,a,n:,,ny,a{t){cr\s\a) , (16) 



where s is the single-particle spin operator. The SRT r 
can then be defined by an exponential fit to the envelope 
of the z component of the average spin of the ensemble 
of electrons)^ 



are taken from the diagonalization of 



Vc{r) and 



(5), = E((s)fe).w- 



(17) 
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In all the numerical results, we include the electron- 
phonon and electron-electron scattering. As initial con- 
dition we assume a spin polarization along the z-direction 
with a small initial spin polarization P — {2 / h) (S) z / 
where Ng is the total electron density. 



TABLE I: Material parameters used in the calculation (from 
Ref . 13 unless otherwise specified) . 



me/ mo 
Vsi (km/s) 
blf' (eA^' 
H(eV) 
Aq (eV) 



12.25 

0.023 

4.28 

27.18" 

5.8 

0.38 



Ko 

Olo (meV) 
Vst (km/s) 
llf ^ (eA) 



ei4 
E, 



(V/m) 
P (eV) 



15.15 

27.0 

1.83 

117.1" 

0.35 X 

0.414 



10" 



" Ref. m 



1200 . . 

= 60 nm — ^ 
= 50 nm ■ 
1000" a, = 40nm . 




FIG. 1: SRT r 
(100) QWRs at different 
4 X 10^ cm"^ and T = 100 K, 



45 50 55 60 65 70 75 80 



ay (nm) 



the QWR width in y direction, 



for 

The electron density is A'' = 
The arrows mark the densities 



at which the electron populations in the second and higher 
subbands are approximately 30 %. 



Influence of the wire size 



In Fig. 1 we plot the SRT as a function of the QWR 
width in y direction, ay, for various a^. We choose (100) 
QWRs at a lattice temperature of T = 100 K and a 
total electron density of Ne = 4 x 10^ cm"-'^. For a 
small QWR of width = 40 nm, the SRT increases 
monotonously with ay because for a small QWR the spac- 
ing of the subbands is large, so that for the present con- 
ditions only the lowest subband in QWR is appreciably 
populated. Therefore, as discussed in Sec. II, the effec- 
tive magnetic field due to the Dresselhaus term contains a 
longitudinal component, B^{k), that keeps the electronic 
spins aligned and thus inhibits spin precession, which is 



mainly due to the Rashba term, thereby effectively re- 
ducing the spin relaxationi^ An interesting effect arises 
because B^{k) is proportional to ((fc^) — {ky)), which 
disappears when the wire widths, and therefore in the 
present model also the confinement wave functions, in x 
and y directions are identical. For vanishing B^{k), i.e., 
for a^ — cLy, the SRT reaches a minimum of several 10 ps 
for the QWRs considered here. Changing the wire size, 
leads to increasing ((fc^) — {ky)) and therefore to increas- 
ing SRT. This is clearly illustrated for the ax = 50 nm 
case. When ax = 60 nm, one can also find a minimum of 
SRT at ax = ay = 60 nm for the reason described above. 
However, one finds that the SRT decreases with ay again 
when ay > 70 nm. This is because the electrons start to 
populate higher subbands when the wire width increases. 
Then the effective magnetic field formed by the Dressel- 
haus term contains not only a longitudinal component 
but also a transverse component, B^{k), which couples 
different subbands. This coupling has the same effect 
as a spin precession, and therefore tends to make the 
SRT shorter. Moreover, the contribution of the Dressesl- 
haus term becomes as important as the contribution of 
the Rashba term. When ay > 70 nm, B^{k) dominates 

over B^ (k) due to the increasing ay, and a faster spin 
relaxation results. The different minima of the SRT at 
o-x — o,y — 60 nm and ax = ay = 50 nm are likely due to 
the different energy gaps between different subbands in 
the two cases. 



B. (110) and (111) QWRs 
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FIG. 2: SRT r vs. for (110) QWRs at T = 50 K and 
Ae = 4 X 10'' cm~^. A: a^, = = 30 nm with an initial spin 
polarization along the a:;-direction; ■: ax = ay = ?,Q nm with 
an initial spin polarization along the y-direction; •: = dy = 
50 nm with an initial spin polarization along the a;-direction. 

The results of the previous subsection showed that 
the geometry of the QWR has a pronounced effect on 
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the SRT because it influences the different contribu- 
tions to the SOC. We now investigate QWRs with dif- 
ferent growth directions, and start with the case of (110) 
QWRs. The SOC for (110) QWR is quite different from 
the SOC for (100) QWR as shown in Eq. ^ and Eq. Q. 
First we only consider the case of a narrow wire with 
dx = ciy = 30 nm, for which the electronic population 
is mainly in the lowest subband. In the presence of an 



electric field of the form (E^ 
butions are 



Ey,0), the relevant contri- 



1 



(18) 
(19) 



The effective magnetic field formed by the Dresselhaus 
term is along the a;-direction, which corresponds to the 
[110] crystallographic direction, and the effective mag- 
netic field formed by the Rashba term is in the x-y 
plane. If the direction of the total effective magnetic field 
formed by the SOC is tuned to be exactly the direction 
of the initial spin polarization, then one can expect an 
extremely long SRT as pointed out in Refs. [H and [H. 
For the QWRs considered in this paper, we can study 
the physics that gives rise to this effect in the following 
way: We take the initial spin polarization to be along the 
x-direction. In Fig. 2 we plot the SRT as a function of 
Ex = EcosB for temperature T = 50 K, electron density 
Nf. = 4x10^ cm~^, and wire geometry ax = cLy = 30 nm. 
The SRT decreases with Ex^ which is a measure of the 
effective magnetic field along y-direction. As expected, a 
very long SRT results when Ex is very small, even though 
the effective magnetic field along the x-direction is not 
zero. In comparison, for an initial spin polarization along 
the y-direction a much shorter SRT is obtained, which 
hardly changes with Ex- This is because the polarized 
electronic spins precess around the x-direction and this 
precession is not influenced strongly by the y-component 
of the effective magnetic fleld. In the inset of Fig. 2, we 
also show the dependence of the SRT on the angle 9, 
which the electric field in the x-y plane forms with the 
X axis. Assuming i? = 10 kV/cm and an initial spin po- 
larization along the x-direction, it is found that the SRT 
decreases with 9. For small 9, i.e., for effective magnetic 
fields close to the x-direction, the SRT goes up, in agree- 
ment with the previous discussion. The result for a larger 
wire size, ax = ay = 50 nm are also plotted in Fig. 2. In 
this case the SRT is never longer than 1 ns, even when 
Ex is very small because now electrons populate higher 
subbands, and the effective magnetic field formed by the 
Dresselhaus term contains not only the longitudinal com- 
ponent but also a transverse component. 

We compare the SRT obtained by our calculation with 
Monte-Carlo results from Ref. [25| in Fig. 3. for a tem- 
perature of T = 5 K. The total electron density and the 
impurity density are taken to be A^e = 4 x 10® cm~^ 
and Ni = 3.17 x 10^ cm~^, and the wire geometry is 
o-x = a,y = 1 fim. The results without the Coulomb 
scattering are actually very close to those obtained by 
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FIG. 3; Our results compare witli the results obtained by 
Monte-Carlo simulations with different wire directions at T = 
5 K, A^e = 4 x 10* cm~\ iV, = 3.17 x lO'' cm"\ and = 
By = 1 ^m. 



Monte-Carlo simulations. However, when the Coulomb 
scattering is included, much longer SRTs result. This re- 
sult underscores the importance of Coulomb scattering 
for spin relaxation in QWRs. 




40 60 80 
(kV/cm) 



100 



FIG. 4: SRT T 



Ex for different wire sizes 



Cby at 



T = 50 K and A^e = 4 x 10^ cm"\ The growth direction of 
the QWR is along the [111] crystallographic direction. The 
solid curves are the results with the initial spin polarization 
along ^-direction and the dashed curve is the result with the 
initial spin polarization along a::-direction. 



For the quantitative analysis of the SRT for another 
growth direction we choose (111) QWRs. Here, we 
again consider the case of small wire width and low tem- 
perature first. The electric field is again taken to be 
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{Ex, Ey,0), so that the SOC can be written as 
Hr^ = -f!f"[-cr,Eyh+ayExk,] , (20) 

-la.k^ikl)}. (21) 

Similar to the case of (110) QWRs, we expect very long 
SRTs if the total effective magnetic field points into the 
direction of the initial spin polarization, at least for nar- 
row wires. For a numerical example of this effect, we 
choose Ey such that j^f^Ey + ( V2/3)62f '^(fc^) = for a 
small wire with = a-y = 10 nm, so that the x compo- 
nent of the total effective magnetic field is zero. For an 
initial spin polarization along the z-direction, which cor- 
responds to the [111] crystallographic direction, we plot 
the SRT as a function of E^ in Fig. 4 at T = 50 K 
and iVe = 4 X 10^ cm~^. Fig. 4 shows that when 
o-x — o,y — 10 nm, there is a pronounced maximum of 
the SRT at E^ = 70 kV/cm, which fulfils the relation 
llf'E, + -l=btf%{kl) + (fc2)) « 0. Consequently, for 

this field strength, the direction of the total effective mag- 
netic field is exactly along the direction of the initial spin 
polarization and this leads to a very long SRT. However, 
this effect much reduced for larger wire cross sections. 
When ax = dy = 20 nm, there is still a maximum of the 
SRT, but the maximum is much less pronounced than 
for the smaller wire, because now "fti^'^Ey + ^h\f"^{ky) 
remains finite for all field strengths as (fc^) is changed. 
This trend continues for wire sizes of ~ ay — 30 nm 
and ax = ay = 50 nm. Finally, we analyze the the case of 
CLx = <iy = 20 nm with the initial spin polarization along 
x-direction: Here one also does not obtain a maximum of 
the SRT because the direction of the effective magnetic 
field is no longer identical to the direction of the initial 
spin polarization. 

C. Doping and temperature dependence 

Since the population of higher subbands has been 
shown to play an extremely important role for the SRT, 
we next analyze how the electronic population can be 
changed without using a different wire geometry or 
growth direction, namely by varying the doping density 
and/or temperature. In Fig. 5(a) we plot the SRT as a 
function of TVg for a (100) QWR of size ax = ay — 50 nm, 
and T = 100 K. For the fuU calculation, the SRT de- 
creases with Nf, because more electrons are present at 
higher momenta and in higher subbands, so that the ef- 
fective magnetic fields experienced by these electrons are 
larger. By the DP mechanism, this results in a faster 
spin relaxationi^ To investigate how scattering affects the 
spin relaxation, we first switch off inter-subband electron- 
phonon scattering, with the result that the SRT becomes 
shorter. When we switch off the inter-subband electron- 
electron scattering, the SRT also becomes shorter. This 
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(a) a,j = By = 50 nm 



full calculation — ■ — 








200 (b) a,^ = ay = 70 nm 

full calculation — ■ — 
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without Coulomb — • — 

without inter-subband electron-phonon — • — 
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5 10 15 20 25 30 
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FIG. 5: SRT r vs. the electron density for (100) QWRs at 
different wire sizes, (a) = Uy = 50 nm; (b) = ay = 
70 nm (b). T = 100 K. 



somewhat counterintuitive effect of scattering on the spin 
dephasing results because the infiuence of scattering is 
different for different regimes, i.e., for strong and weak 
scattering: If we define — Hr + Ho and r* to be the 
effective momentum relaxation time, then |ri|r* is typ- 
ically much smaller than 1 for the densities considered 
here {\n\T* = 0.03 at iVg = 10 x 10^ cm^^ if t* contains 
only the contribution from the inter-subband electron- 
electron scattering). This means we are in the strong 
scattering regime where the motional narrowing picture 
qualitatively describes the dependence of the SRT on the 
momentum scattering time as SRT w 1 /r* . (Microscop- 
ically this means that the scattering is strong enough 
to prevent significant deviations from isotropic electron 
distributions, which are caused by the anisotropic SOC 
contributions.) Thus switching off scattering contribu- 
tions effectively leads to a longer momentum relaxation 
time, and therefore to shorter SRTs, while still remaining 
in the strong scattering regime. Further, our calculation 
shows that the inter-subband electron-electron scatter- 
ing is more important than the inter-subband electron- 
phonon scattering. When we further exclude the intra- 



7 



subband electron-electron scattering, it is found that the 
SRT becomes shorter when A^e < 20 x 10^ cm"-^, but 
becomes longer when A'e > 20 x 10^ cm^^. Here the 
most dominant scattering process is suppressed, so that 
with increasing Nf, the motional narrowing regime is left, 
because electrons occupy states at higher momenta and 
higher subbands where the SOC contributions rapidly in- 
crease. For higher densities, the weak scattering regime 
is reached where the anisotropy of the SOC contributions 
becomes dominant. Additional scattering leads to a more 
efficient dephasing, so that the SRT increases when we 
switch off scattering contributions. 

We also plot the case of larger wire width with = 
Qy = 70 nm in Fig. 5(b). Compared to the case of 
Fig. 5(a), the electrons are populating higher subbands, 
so that the Dresselhaus term becomes more important. 
Therefore, the band structure anisotropy will become 
important at smaller Ne- Fig. 5(b) shows that the re- 
sult of the full calculation intersects the one without 
the electron-electron scattering at A^g « 8 x 10^ cm~^. 
For this wire geometry, the result without the inter- 
subband electron-electron scattering also crosses the full 
calculation. To see in more detail how the strength of 
anisotropic SOC contributions affects the spin relaxation, 
we also plot the result without the Dresselhaus term. 
For this somewhat artificial case, the SRT becomes much 
longer. This indicates that the contribution of the Dress- 
eslhaus term is very important for this wire width in con- 
trast to the case of smaller wire width in Fig. 5(a), be- 
cause the transverse component of the effective magnetic 
field contributed by the Dresselhaus term is important 
only for the large wire width when higher subbands are 
populated. 
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FIG. 6: SRT T vs. the temperature for (100) QWRs at differ- 
ent wire sizes ax ~ ay at carrier density A*" = 4 x 10^ cm~^. 

Finally we look at the temperature dependence of 
the SRT at different wire widths in Fig. 6. We choose 
(100) QWRs at a representative electron density of A^^ = 
4 X 10^ cm~^. The temperature affects the SRT in two 



ways: For the smaller wires = ciy < 50 nm the carri- 
ers are confined in the lowest subband without populat- 
ing higher fc-states where the anisotropic SOC contribu- 
tions to the bandstructure become much stronger. In this 
strong scattering regime, the behavior can qualitatively 
be explained by the motional narrowing: Increasing tem- 
perature leads to enhanced scattering with phonons and 
electrons with higher kinetic energies, so that the effec- 
tive carrier lifetime becomes shorter, and the SRT in- 
creases. Around ax = ay = 60 nm and T > 150 K 
this behavior crosses over to the weak scattering limit, 
in which higher subbands are populated, the precession 
frequencies around the effective internal fields become 
higher, and additional scattering leads to a more efficient 
dephasing, so that the SRT decreases with T. 



IV. CONCLUSION 

In conclusion, we have investigated the spin relaxation 
of electrons in n-type InAs QWRs. The SRT is calcu- 
lated by numerically solving the microscopic KSBEs in- 
cluding multiple subbands. The inclusion of higher sub- 
bands allows us to investigate QWRs larger QWRs than 
in Ref. ,3L . and we find that the quantum- wire size influ- 
ences the spin relaxation time via the SOC: The Dres- 
selhaus term contains a longitudinal contribution to the 
internal effective magnetic field, which can effectively re- 
duce the spin precession, and thereby the spin relaxation. 
It also contains a transverse component, which provides 
a contribution to the spin precession involving different 
subbands and enhances spin relaxation. When the wire 
width is small and only the lowest subband in the QWR 
is important, the longitudinal term is dominant. When 
the wire width is large and higher subbands are pop- 
ulated, the transverse contribution dominates over the 
longitudinal one, and a faster spin relaxation results. We 
also studied different growth directions for QWRs. We 
show that one can obtain long spin relaxation time by op- 
timizing the growth direction, quantum-wire width and 
the direction of the initial spin polariazation. Further, we 
investigated how the details of the microscopic scattering 
mechanisms and the spin-orbit effects in the band struc- 
ture affect the spin relaxation. The population of higher 
subbands was found to have decisive influence on the be- 
havior of the SRT. For instance, if the geometry and ex- 
ternal conditions are such that higher subbands become 
populated, the dependence of the SRT on temperature is 
reversed because the motional narrowing regime is left. 
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